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Whole-genome duplication events have driven to a large degree the evolution of angiosperm genomes. 
Although the majority of redundant gene copies after a genome duplication are lost, subfunctionalization or 
gene balance account for the retention of gene copies. The Arabidopsis SOS ribosome represents an excellent 
model to test the gene balance hypothesis as it consists of 80 ribosomal proteins, all of them encoded by 
genes belonging to small gene families. Here, we present the isolation of mutant alleles of the APICULATA2 
{API2) and RPL36aA paralogous genes, which encode identical ribosomal proteins but share a similarity of 
89% in their coding sequences. RPL36aA was found expressed at a higher level than API2 in the wild type. 
The loss-of- function api2 and rpl36aa mutations are recessive and affect leaf development in a similar way. 
Their double mutant combinations with asymmetric leaves2-l {as2-l) caused leaf polarity defects that were 
stronger in rpl36aa as2-l than in api2 as2-L Our results highlight the role of combined haploinsufficiency 
and purifying selection in the retention of these paralogous genes in the Arabidopsis genome. 



Ancient genome duplication (polyploidization) events have played a central role in the evolution of 
angiosperm genomes^ Several models account for the fate of individual genes after a duplication occurs. 
The redundant gene copies might be inactivated and lost or, alternatively, they might be retained in the 
genome by different evolutionary processes. In the duplication -degeneration -complementation model, subfunc- 
tionalization (i.e. the distribution of the ancestral functions among the duplicated genes) accounts for the 
evolutionary maintenance of duplicate gene copies^. Duplicated genes might also be maintained, retaining the 
same (redundant) function and expression pattern, when the duplication leads to a new balance between 
the products of dosage-dependent genes (gene balance hypothesis) 

The availability of complete plant genome sequences has revealed the presence of recent and ancient whole- 
genome duplication events in the evolutionary history of numerous plant species. Three rounds of polyploidiza- 
tion, followed by extensive chromosomal rearrangements and gene loss to reach the present diploid state, have 
been documented in the genome of Arabidopsis thaliana (hereafter, Arabidopsis) ^ The phenomenon by which 
individual genes are lost from one of the two homeologous chromosome pairs after a whole-genome duplication 
event is called fractionation^. Somewhat counter- intuitively, it has been demonstrated that fractionation is a 
biased phenomenon, with a majority of genes being preferentially lost from one of the homeologous chromo- 
somes^. Interestingly, some genes appear to have escaped fractionation, so that two or more gene copies resulting 
from duplication events have persisted in the genome^. While retention of some gene pairs might be attributed to 
subfunctionalization, the gene balance hypothesis accounts for the retention in other cases. Indeed, the retained 
genes have been found to preferentially belong to certain over-represented functional categories^, which typically 
comprise proteins that are dosage-dependent or form part of multisubunit complexes, often with regulatory roles, 
suggesting that haploinsufficiency has played a major role in their retention. Loss-of-function mutations at 
haploinsufficient loci are typically dominant because the level of gene function in a heterozygote is below the 
threshold to produce a wild-type phenotype. 

The eukaryotic ribosome represents an excellent model for testing the gene balance hypothesis, as it consists of 
a large subunit (60S) with 46 ribosomal proteins and 3 ribosomal RNA (rRNA) molecules (25S, 5.8S and 5S), and 
a small subunit (40S) with 33 ribosomal proteins and a single 18S rRNA molecule^. The structure of the 
prokaryotic and eukaryotic ribosomes has been determined at an atomic-level resolution^^"^^. Prokaryotic and 
eukaryotic ribosomes share a general macroorganization but they differ in size and number of ribosomal proteins 
and rRNA molecules. In Arabidopsis, all cytosolic ribosomal proteins are encoded by small gene families, each 
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with 2 to 7 members^^. Although 249 genes encoding ribosomal 
proteins have been annotated in the Arabidopsis genome, a max- 
imum of 80 different ribosomal proteins has been identified in iso- 
lated SOS ribosomes^^"^^. 

Here, we report the isolation of mutant alleles of two genes, 
APICULATA2 {API2) and its closest paralog (i.e. homologous gene 
resulting from a duplication event) in the Arabidopsis genome, 
At3g23390. These genes encode two identical ribosomal proteins, 
RPL36aB and RPL36aA, despite they carry nucleotide substitutions 
at 34% of their codons. Loss-of-function mutations at each of these 
unlinked loci are recessive and affect leaf development in a similar 
way. We found that the api2 and rpl36aa mutations represent a rare 
example of non-allelic non-complementation, indicating that two 
(out of the four) functional gene copies present in a diploid plant 
are not sufficient to produce a wild-type phenotype. We demonstrate 
that both paralogous genes are expressed in the same tissues and 
show that both loci behave as a single haploinsufficient unit with 
four alleles (combined haploinsufficiency) on which natural selec- 
tion has acted to filter out deleterious non- synonymous substitutions 
(i.e. nucleotide substitutions that change the protein sequence), a 
process known as purif)^ing selection. 

Results 

Cloning and molecular characterization of the api2 mutation. In a 

large-scale screen for ethylmethane sulfonate (EMS) -induced mu- 
tants with abnormal leaf shape, we previously isolated 42 mutants 
with a pointed-leaf phenotype. These mutants were assigned to 
phenotypic classes dubbed Apiculata (Api), Denticulata (Den) and 
Angusta (Ang)^°'^\ The Api class comprises recessive alleles of seven 
different genes, named APIl to API7^^. Here, we focus on the API2 
gene, which was initially defined by a single mutant allele, api2. To 
identif)^ API2 at the molecular level, we undertook a positional 
approach. We first mapped the api2 mutation to an 850-kb in- 
terval on chromosome 4, flanked by two insertion-deletion (InDel) 
markers, cer459267 and cer453616 (Figure la), encompassing 248 
annotated genes. Because the api2/api2 mutant has pointed leaves 
(Figure 2b, d, g, i), and is very similar in phenotype to other 
mutants affected in genes that encode ribosomal proteins^^"^^, we 
focused on the only two genes within the interval that encode 
subunits of the cytosolic ribosome, At4gl3170 and At4gl4320. After 
sequencing the transcriptional unit of both genes, we found a single G- 
to-A transition mutation in the 5' untranslated region (UTR) of the 
At4gl4320 gene (Figure lb and Supplementary Figure SI), also 
known as RPL36aB^^. This mutation creates an out-of-frame ATG 
codon upstream of the normal translation start site of RPL36aB 
transcripts (Supplementary Figure SI). We used real-time qRT-PCR 
to measure the relative RPL36aB expression level in rosettes collected 
14 days after stratification (das), and found a 0.57-fold decrease in 
api2/api2 plants compared with the wild type (Figure 3). To deter- 
mine whether the lesion in RPL36aB is responsible for the mutant 
phenotype, we followed a transgenic approach to complement the 
api2 mutation. The phenotype of api2/api2 plants expressing a 
35Spro:RPL36aB transgene was indistinguishable from wild type in 4 
independent lines (Figure 21, o), indicating that loss of RPL36aB 
ftmction causes the mutant phenotype, either by the observed 
down -regulation of transcript levels or because the api2 mutation 
interferes with normal translation initiation. 

Two paralogous genes encode the RPL36a protein in Arabidopsis. 

The Arabidopsis genome contains two RPL36a paralogous genes that 
encode identical proteins, RPL36aA (At3g23390) and the above- 
mentioned RPL36aB (hereafter, API2y^. We also searched for 
mutant alleles of RPL36aA, and identified a publicly available T- 
DNA insertion line (SALK_148438) that, when homozygous for 
the insertion, exhibited a pointed-leaf phenotype identical to that 
caused by the api2 mutation in the Col-0 genetic background 
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Figure 1 | Positional cloning of the API2 gene and molecular nature and 
position of the api2 and rpl36aa mutations, (a) A mapping population of 
138 F2 plants derived from an api2/api2 X Col-0 outcross was used to 
define a candidate interval of 850 kb on chromosome 4. The names and 
map positions of the molecular markers used are indicated. The number of 
recombinant chromosomes found and the total of chromosomes analyzed 
are shown in parentheses. Black boxes represent BAG clones, (b) Structure 
of the AP12 and RPL36aA paralogous genes, with indication of the nature 
and position of the api2 and rpl36aa mutations. A triangle indicates the 
position of a T-DNA insertion. Boxes and lines between boxes indicate 
exons and introns, respectively. Open boxes correspond to the 5' and 3' 
untranslated regions. The predicted translation start (ATG) and stop 
(TAA) codons are indicated. 

(Figure 2d, e, i, j). We refer to this line as rpl36aa/rpl36aa. The 
presence of the T-DNA insertion in the first intron of RPL36aA 
was confirmed by PGR genotyping (Figure lb). The inflorescences 
of rpl36aa/rpl36aa and api2/api2 plants were shorter than those of 
wild-type plants (Figure 2k). We used real-time qRT-PGR to 
measure the relative RPL36aA expression level in rpl36aa/rpl36aa 
rosettes collected 14 das, and only found a residual expression level 
(Figure 3), suggesting that rpl36aa is a null allele. To demonstrate 
that the insertion causes the defects observed in the rpl36aa line, 
we also followed a transgenic approach to complement its mutant 
phenotype. Four independent rpl36aa/rpl36aa;35Spro:RPL36aA 
transgenic lines exhibited a wild- type leaf phenotype (Figure 2m, 
p), confirming that the mutant phenotype results from the inactiva- 
tion of RPL36aA. We concluded that identical recessive leaf mutant 
phenotypes occur when either API2 or RPL36aA are inactivated. 

API2 and RPL36aA are dosage-dependent genes. Because the API2 
and RPL36aA proteins have identical amino acid sequences 
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Figure 2 | Phenotypes and phenotypic rescue of the api2 and rpl36aa mutants, (a-e), (1-n) Rosettes and (f-j), (o-q) first-node leaves from (a), (f) Ler, 
(b), (g), (d), (i) apil/apil (b), (g) in a Ler background and (d), (i) in a Col-0 background, (c), (h) Col-0, (e), (j) rpl36aal rpl36aa, (1), (o) AP12/ 
api2;RPL36aA/rpl36aay (m), (p) api2/ api2;35Spro:RPL36aB and (n), (q) rpl36aa/ rpl36aa;35Spro:RPL36aA plants. Arrowheads indicate marginal teeth in 
mutant pointed leaves, (k) Adult plants. All plants are in a Col-0 background, the only exceptions being (b), (g) the api2 mutant and (a), (f) the Ler wild 
type. Pictures were taken (a-j), (1-q) 14 and (k) 35 das. Scale bars indicate (a-j), (1-q) 2 mm and (k) 3 cm. 



(Supplementary Figure S2), they are expected to perform identical 
functions in the ribosome, offering a unique opportunity to ask 
questions about the functional redundancy and haploinsufficiency 
of paralogous ribosomal proteins, with no interference due to the 
divergence of their sequences. The api2/api2 and rpl36aa/rpl36aa 
single mutants are viable, suggesting that both genes contribute 
quantitatively to a global pool of RPL36a ribosome subunits and 
that two functional copies of any one gene can alleviate the effects 
of having two inactive copies of its paralog. We hypothesized that the 
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Figure 3 | Expression of the API2 and RPL36aA genes in the Col-0 wild 
type and the api2/api2 and rpl36aa/rpl36aa mutants. Real-time qRT-PCR 
analysis of the expression of the API2 and RPL36aA genes was performed in 
api2/api2 and rpl36aa/rpl36aa vegetative leaves. Relative expression levels 
are shown, determined by the comparative Cj method, and normalized 
with the expression of the 18S rRNA housekeeping gene. Error bars 
indicate the interval delimited by 2~^^^^x ~ Asterisks indicate ACj 
values significantly different from those of Col-0 in a Mann-Whitney U test 
(p < 0.01; n = 9). 



simultaneous loss of one copy of each locus should also lead to a 
mutant phenotype similar to those of api2/api2 or rpl36aa/rpl36aa 
single mutants. We crossed the single mutants and found that all the 
doubly heterozygous plants (API2/api2;RPL36aA/rpl36aa) in the Fi 
progeny exhibited a pointed-leaf phenotype identical to those of the 
api2/api2 and rpl36aa/rpl36aa homozygous parentals (Figure 2n, q). 
This non- allelic non- complementation involving hypomorphic^^ 
(partial loss-of-function) alleles of API2 and RPL36aA requires 
that both genes play redundant functions and are co-expressed in 
at least some cells. 

In order to investigate the effects of carrying less than two wild- 
type alleles, we genotyped 169 mutant plants from an F2 family 
segregating mutant alleles of both loci. We did not find plants with 
a single functional copy of either gene {API2/api2;rpl36aa/rpl36aa or 
api2/api2;RPL36aA/rpl36aa) or lacking wild-type copies {api2/ 
api2;rpl36aa/rpl36aa), despite these genotypes are collectively 
expected to appear in 3L25% (5 out of 16) of the F2 plants. 
Interestingly, the progeny of self-pollinated API2/api2;RPL36aA/ 
rpl36aa plants segregated an excess of aborted ovules (Figure 4), with 
a 13% increase relative to the single mutants and the wild type. 

To test any male- or female-gametophytic lethality, we per- 
formed reciprocal crosses between doubly heterozygous API2/api2; 
RPL36aA/rpl36aa and Col-0 plants. When we pollinated API2/ 
api2;RPL36aA/rpl36aa pistils with Col-0 wild-type pollen, we 
observed an average 23% excess of aborted ovules in the siliques 
relative to the values obtained in control crosses, indicating female 
gametophytic lethality. Reciprocally, we also pollinated Col-0 pistils 
with pollen from API2/api2;RPL36aA/rpl36aa plants. If the four 
possible types of haploid male gametophytes are equally transmitted 
to the progeny, then one fourth of the Fi plants should have the same 
genotype {API2/api2;RPL36aA/rpl36aa) and exhibit the same 
mutant phenotype as their male parental. However, all 140 Fi plants 
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Figure 4 | Aborted ovules in API2/api2;RPL36aA/rpl36aa siliques. 

(a-d) Dissected mature siliques from (a) Col-0, (b) apillapil (Col-0), 
(c) rpl36aa/rpl36aa and (d) AP12I api2;RPL36aAI rpl36aa plants, 
(e), (f) Detail of an aborted ovule in an API2I api2;RPL36aAI rpl36aa 
cleared silique. Arrowheads in (d) and (e) indicate aborted ovules. Scale 
bars indicate (a-d) 1 mm and (e), (f) 100 |am. 

studied showed a wild-type phenotype, a result that indicates that the 
api2 rpl36aa gametophytes are not viable. 

RPL36a genes are subjected to strong purifying selection. The 

identical protein sequences encoded by different RPL36a paralo- 
gous genes might occur as a consequence of their recent evolu- 
tionary origin; alternatively, they might be explained by concerted 
evolution (e.g. gene conversion) or by strong purifying selection^^. 
The coding regions of both API2 and RPL36aA transcripts consist of 
105 codons. A pairwise alignment of the API2 and RPL36aA coding 
sequences uncovered 36 nucleotide substitutions, each affecting a 
different codon (Supplementary Figure SI). All these substitutions 
were synonymous, affecting the third base of a codon in 35 cases and 



the first base only in 1 case (Supplementary Figure S2). Because the 
affected codons were scattered along the sequence, this level of 
divergence is unlikely to have arisen as a consequence of concerted 
evolution or a recent duplication. We determined the rates of 
synonymous (ds) and non-synonymous substitutions (d^) from a 
pairwise alignment of the API2 and RPL36aA coding sequences. The 
low value of the d^/ds ratio (Table 1) is a signature of strong 
purifying selection, and is in line with the values obtained for other 
pairs of paralogous genes encoding identical ribosomal proteins. 

API2 and RPL36aA are expressed in overlapping patterns. As 

mentioned above, the non-allelic non- complementation between 
API2 and RPL36aA is best understood when there is at least some 
overlap of expression in the same cells. To experimentally demon- 
strate this, we fused the promoter regions of both genes to the GUS 
reporter gene. GUS staining of 5 APHp^o'-GUS and 5 RPL36aAp^o- 
GUS independent transgenic lines in the Col-0 background 
uncovered an identical expression pattern for both genes. Both 
constructs were ubiquitously expressed, with particularly intense 
GUS staining in young and actively proliferating tissues, such as 
those of developing leaves, floral buds and root apices (Figure 5). 
Interestingly, the relative intensities of the GUS staining matched our 
qRT-PCR results (Figure 3), with the RPL36aA promoter conferring 
higher expression levels in wild-type vegetative leaves than the API2 
promoter. 

api2 and rpl36aa genetically interact with asymmetric leaves 
mutants. Some ribosomal proteins are known to be involved in 
leaf developmental programs, as their mutations enhance the 
adaxial-abaxial polarity defects of the asymmetric leaves 1 (asl) and 
as2 mutants^'*'^^ Moreover, Horiguchi et al?^ showed that the extent 
of this enhancement varies among mutants affected in different 
ribosomal proteins. To determine whether api2 and rpl36aa can 
enhance the polarity defects of as mutants, we isolated the api2 
asl-1, api2 as2-l, rpl36aa asl-1 and rpl36aa as2-l double mutants. 
Because both paralogous proteins share an identical sequence and 
expression pattern, any phenotypic differences observed in these 
mutants are likely to reflect the contribution of individual paralogs 
to the pool of RPL36a subunits. Both api2 and rpl36aa enhanced the 
phenotype of the as2-l mutant, although to a different extent 
(Figure 6e-h). Both api2 as2-l and rpl36aa as2-l double mutants 
exhibited trumpet- shaped (peltate) leaves (Figure 6f, g), which we 
did not observe in the as2-l single mutant (Figure 6d, e). In addition 
to this, the rpl36aa as2-l double mutant exhibited radial leaves 
(Figure 6h), which can be interpreted as a consequence of a more 
severe abaxialization. These phenotypes were not observed in 
the double mutant combinations involving the asl-1 mutation 
(Figure 6b, c). 



Table 1 | Selected ribosomal protein paralogs and their synonymous and non-synonymous substitution rates 







Protein length 


Amino acid 


Nucleotide sequence 








Para log pairs 


AGI gene codes 


(aa)* 


sequence identity (%) 


identity (%) 


ds 


dn 


dn/ds 


API2, RPL36aA 


At4g 14320, At3g23390 


105 


100 


89 


1.1338 


0.0011 


0.0010 


RPS18A, RPS18C 


Atlg22780, At4g09800 


152 


100 


88 


1.1036 


0.0056 


0.0051 


RPS28A, RPS28B 


At3g 10090, At5g03850 


64 


100 


88 


0.9059 


0.0009 


0.0010 


RPL21A, RPL21D 


Atlg57860, Atlg57660 


164 


100 


99 


0.0339 


0.0000 


0.0010 


RPL38A, RPL38B 


At2g43460, At3g59540 


69 


100 


90 


0.7566 


0.0008 


0.0010 


RP141A, RPL41B 


At2g40205, At3g08520 


25 


100 


95 


0.3252 


0.0003 


0.0010 


RPS13ARPS13B 


At4g00100, At3g60770 


151 


99 


85 


1 .5649 


0.0029 


0.0018 


RPS21A, RPS21B 


At3g53890, At5g27700 


82 


94 


83 


3.9235 


0.0297 


0.0076 


RPlUA, RPlUC 


At2g37190, At5g60670 


166 


93 


81 


1.01 14 


0.0572 


0.0566 


RPL28A, RPL28B 


At2g 19730, At4g29410 


143 


90 


85 


0.7896 


0.0597 


0.0756 



*AII paralog pairs shown in this table have identical protein length. 
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Figure 5 | Visualization of API2pj-o:GUS and RPL36aApj.o:GUS transgene activity in the Col-0 background, (a), (m) Seedlings and (b), 

(n) cotyledons, (c-f), (o-r) Vegetative leaves from the (c), (o) first, (d), (p) third (e), (q) fifth and (f), (r) seventh nodes, (g), (s) Inflorescences. 

(h), (t) Immature and (i), (u) mature flowers, (j), (v) Immature siliques, (k), (w) root tips and (1), (x) mature siliques. Pictures were taken (a), (m) 6, 

(b-f), (k), (n-r), (w) 14 and (g-j), (1), (s-v), (x) 35 das. Scale bars indicate (a), (c-e), (g), (j), (1), (m), (o-q), (s), (v), (x) 1 mm, (b), (f), (i), (n), (r), (u) 

500 i^m, (h), (t) 200 [im and (k), (w) 100 |^m. 



Discussion 

Following a positional cloning approach, we have determined that 
the api2 mutant carries a lesion at a gene encoding the RPL36a 
ribosomal protein. This ribosomal protein is encoded by two para- 
logous genes, API2 {RPL36aB) and At3g23390 {RPL36aA), which 
have not previously been characterized in Arabidopsis. Our results 
indicate that both paralogs are functionally equivalent (redundant), 
as shown by the non-allelic non-complementation observed in the Fi 
progeny of crosses involving mutant alleles of both genes. This func- 
tional equivalence is also supported by the identical protein se- 
quences and expression patterns conferred by the promoter regions 
of API2 and its paralog. This non-allelic non- complementation indi- 
cates that a critical level of protein function is not reached in doubly 
heterozygous plants, a quantitative effect that is expected only when 
the two genes are expressed in the same cells and perform identical 
function. Although this type of non-complementation is a relatively 
rare phenomenon, numerous additional instances have been docu- 
mented in Arabidopsis, often involving ribosomal proteins of the 
large and small subunits of the cytosolic ribosome^^'^^'^°'^\ The obser- 
vation that loss-of-function mutations damaging ribosomal proteins 
are usually recessive in Arabidopsis^^, with only a few examples of 
semi-dominance^^'^^, is in striking contrast with the dominance or 
haplolethality of loss-of-function (haploinsufficient) Minute muta- 
tions of Drosophila melanogaster^'^ . The Minute mutations represent 
the most abundant phenotypic class in Drosophila and affect genes 
coding for ribosomal proteins. However, unlike in Arabidopsis, each 
ribosomal protein is typically encoded by a single gene in the genome 
of Drosophila^^. This fundamental difference reflects how the need to 
keep a balance between the different ribosomal proteins after a 
whole-genome duplication has shaped the genome during plant 
evolution. 



Loss of API2 or RPL36aA functions led to a deleterious recessive 
syndrome that was very similar, not to say undistinguishable, from 
the phenotypes of many other mutants carrying lesions in genes 
encoding ribosomal proteins^^'^^"^^'^^'^^'^^. Our reciprocal crosses indi- 
cate that the simultaneous inactivation of both paralogs causes leth- 
ality and prevents transmission of the mutant alleles through the 
male and female gametophytes, as has also been reported for other 
pairs of paralogous ribosomal proteins in Arabidopsis^^. Indeed, we 
were not able to isolate plants with fewer than two active (out of the 
four) copies of the genes encoding the RPL36a ribosomal protein. 
The observation that only synonymous nucleotide substitutions have 
accumulated in the coding region of these paralogous genes (none of 
which changes the protein sequence) suggests that a matching num- 
ber of non -synonymous changes have also occurred but must have 
been selected against by strong purifying selection. A similar situ- 
ation is likely to hold for other ribosomal proteins for which lethality, 
non-allehc non- complementation and/or dosage effects have been 
described. The presence of identical sequences in both paralogs (at 
the amino acid level) is not a peculiarity of the RPL36a paralogy 
group, as several instances of gene pairs encoding identical ribosomal 
proteins have been reported in the literature (e.g. RPS18, RPS29, 
RPS30, RPLll, RPL21, RPL23, RPL38, RPL40 and RPL41)'l We 
have found that these proteins have also been subjected to strong 
purifying selection, as inferred from the low d^/ds values for pairs of 
paralogous genes, highlighting the critical roles played by such pro- 
teins in the ribosome. Purifying selection has also been shown to 
account for the sequence conservation of other proteins that play 
essential, highly conserved functions in eukaryotes, such as members 
of the histone H4 protein family^^. We propose that combined hap- 
loinsufficiency and strong purifying selection have driven the biased 
retention of genes coding for paralogous ribosomal proteins after the 
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Figure 6 | Double mutant combinations of api2 and rpl36aa with asl-1 
and as2-l. (a-d) Rosettes of the (a) asl-1, (b) apz2 asl-l, (c) rpl36aa asl-1 
and (d) as2-i mutants, (e-f) Details of vegetative leaves from (e) as2-l, 
(f) api2 as2-l and (g), (h) rpl36aa as2-l mutant plants. All these mutants 
were homozygous for the mutations indicated. Pictures were taken (a-c) 
14 and (d-f) 18 das. Scale bars indicate 1 mm. 

duplication of plant genomes, as has been documented for the 
Arabidopsis genome^ ^ The stoichiometry and structural constrains 
of a multisubunit complex, such as the ribosome, are likely the ulti- 
mate causes behind this retention^^. 

Because API2 and RPL36aA are identical at the amino acid level, 
any phenotypic difference between the api2 and rpl36aa loss-of- 
function mutant alleles should be attributed to differences in the 
spatio-temporal expression pattern or expression levels of the two 
genes. In our particular case, we observed that rpl36aa enhanced the 
polarity defects of the as2-l mutant (as inferred from the fully radia- 
lized lateral organs seen in the rpl36aa as2-l double mutant) to a 
greater extent than api2 (as the api2 as2-l double mutant showed 
trumpet-shaped, but not fully radialized leaves). Interestingly, we 
found RPL36aA to be expressed at higher levels than API2 in wild- 
type plants, suggesting that the stronger effect caused by the rpl36aa 
mutation is due to the fact that RPL36aA is the main contributor of 
RPL36a subunits. Our results favor the hypothesis that, at least in 
plants, the phenotypic differences observed between mutant alleles of 
paralogous genes, which have been reported in some cases^^'^^, might 
result from differences in the expression patterns or the expression 



levels of individual paralogous proteins, rather than from the differ- 
ential specificity for certain mRNAs of an heterogeneous population 

of ribosomes^^'^^'^^"^\ 

Methods 

Plant material and growth conditions. Arabidopsis thaliana (L.) Heynh. 
Columbia-0 (Col-0) and Landsberg erecta (Lev) wild-type accessions and the 
SALK_148438 line (named here rpl36aa; accession number N648438), were obtained 
from the Nottingham Arabidopsis Stock Centre (NASC). The api2/api2 mutant was 
isolated after EMS mutagenesis of Landsberg erecta (Lev) seeds^°. To standardize 
genetic backgrounds, we outcrossed the api2/api2 mutant three times to Col-0 [api2/ 
api2 (Col-0)]. Growth conditions and crosses were performed as described 
previously^"'^^. 

Map-based cloning. The api2 mutation was low- resolution mapped as described in 
Ponce et aV^. For fine-mapping, we used the cer466314, cer459267 and cer453616 
InDel polymorphisms (Table SI) from Monsanto (http://www.arabidopsis.org/ 
browse/Cereon). To identif)^ the api2 mutation, the entire transcriptional unit of 
At4gl4320 was amplified using the API2_F1 and API2_R1 primers (Table SI). The 
resulting 1660-bp PCR product was sequenced on an ABI PRISM 3130x1 Genetic 
Analyser (Applied Biosystems). The presence of the T-DNA insertion in rpl36aa was 
confirmed by PCR using the rpl36aa_F/R primers designed with the T-DNA Primer 
Design tool (http://signal.salk.edu/tdnaprimers.2.html) (Table SI). 

RNA isolation, cDNA synthesis and qRT-PCR. Total RNA from 20-30 mg of 
rosette leaves, collected 14 das, was isolated using TRI Reagent (Sigma) and first- 
strand cDNA synthesis and relative expression quantification were performed as 
described in Jover-Gil et al.'^'^, using the primers listed on Table SI. 

Genetic constructs and plant transformation. To make the 35Spro:RPL36aB and 
35Spro:RPL36aA transgenes, we amplified the full-length coding sequence of 
At4gl4320 and At3g23390 genes from Col-0 cDNA using Phusion polymerase 
(Thermo Scientific) and the primers shown in Table SI. The resulting products were 
cloned into the pENTR/D-TOPO vector (Invitrogen) and transferred into the 
pMDC32 destination vector^^ in a reaction with LR clonase II (Invitrogen). To make 
the RPL36aBpro:GUS and RPL36aApro.GUS transgenes, 389-bp and 716-bp fragments 
from the upstream region of the starting codon of At4gl4320 and At3g23390 genes, 
respectively, were amplified from Col-0 genomic DNA using Phusion polymerase 
and the API2pro_F/R and RPL36aApro_F/R primers, which include the attBl and 
attB2 sequences at their 5' ends (Table SI). The amplification products were cloned 
into the pGEM-T Easy221 entry vector (kindly provided by B. Scheres), by 
recombination with BP clonase II (Invitrogen), and then transferred into pMDC163 
destination vector"^^ in an LR clonase II (Invitrogen) reaction. Competent 
Agrobacterium tumefaciens LBA4404 cells were transformed with these constructs. 
Col-0, api2 and rpl36aa plants were transformed by the floral dip method^*^. 

GUS assay and microscopy. GUS staining was performed as described in Robles et 
aV, with minor modifications. In brief, plant tissue was incubated in ice-cold 90% 
acetone for 15 min and then in X-Gluc buffer solution (2 mM 5-bromo-4-chloro-3- 
indolyl-P-glucuronic acid, 50 mM sodium phosphate, pH 7.2, 5 mM potassium 
ferrocyanide, 50 mM potassium ferricyanide, and 0.2% Triton X-100) for 14 h at 
room temperature. After GUS detection, the tissue was cleared with 96% ethanol and 
mounted in an 8 : 2 : 1 (chloral hydrate : glycerol : water) solution. Samples were 
examined with a Nikon D-Eclipse CI microscope. 

Bioinformatic analysis of selection. To determine the impact of selection, we 
calculated synonymous (ds) and non- synonymous (dN) substitution rates for pairs of 
aligned sequences using the codeml program of the PAML (Phylogenetic Analysis by 
Maximum Likelihood) package^**, as implemented in the PAL2NAL program (http:// 
www.bork.embl.de/pal2nal/)'^''. 
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